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Abstract. The relationship between the compressible magnetohydrodynamic flows 
with low Mach number and the incompressible magnetohydrodynamic flows is investi- 
gated. More precisely, the convergence of weak solutions of the compressible isentropic 
viscous magnetohydrodynamic equations to the weak solutions of the incompressible 
viscous magnetohydrodynamic equations is proved as the density becomes constant 
and the Mach number goes to zero, that is, the corresponding incompressible limits are 
justifled when the spatial domain is a periodic domain, the whole space, or a bounded 
domain. 



1. Introduction 

Studies on magnetohydrodynamic flows always involve a choice at the onset to de- 
scribe the system entirely in the context of either incompressible magnetohydrodynamics 
(MHD), or compressible MHD. For example, theoretic studies on turbulence have a partic- 
ular leaning toward the incompressible model. This preference has largely been based on 
the benefits and advantages of the similarity of incompressible MHD to its hydrodynamic 
counterparts, and the practical consideration of limited computational resources. However, 
when the density of a flow is no longer invariant, the flow become much more complicated 
not only from the physical viewpoint, but also from the mathematical consideration, see 
131 [IHl [IZl [m [22] and references therein. Thus, it is a natural problem to consider the 
relation between the incompressible MHD and the compressible MHD. The equations of 
the isentropic compressible viscous magnetohydrodynamic flows in N spatial dimensions 
have the following form ( [1 [21] [2^ ) : 

Pt + div(pu) =0, 

{pu)t + div {pu (K) u) -f Vp(p) = (V X H) X H + ptAu + AVdivu, (1.1) 
Ht - V X (S X H) = -V X (j;v X H), divH = 0, 

where p. > is the shear viscosity, A is the bulk viscosity satisfying 2/1 + A > , u > 
is the magnetic viscosity; and p denotes the density, u G the velocity, H € the 
magnetic field, p{p) = ap'^ the pressure with constant a > and the adiabatic exponent 
7 > 1. The symbol (g) denotes the Kronecker tensor product. The first equation in p.ip 
is called the continuity equation and the third equation in (|l.ip is called the induction 
equation. 

From the physics point of view, the compressible flow behaves asymptotically like an 
incompressible flow when the density is almost constant, and the velocity and the magnetic 
flcld are small, in a large time scale. More precisely, we scale p, u, and H in the following 
way: 

p=:p{x,et), u = eu{x,et), ll — eil{x,et), (1.2) 
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and we assume that the coefficients /i, A, and v are small and scaled as: 

Jl = ep,£, X = eXe, V ~ eve, (1-3) 

where e G (0; 1) is a small parameter and the normalized coefficients p,^, Xg and satisfy 

^ /i, X^ ^ A, ^ V, as £ ^ 0+, (1-4) 

with jjL > 0, 2/i + A^A > 0, andj^ > 0. Such a scaling as (|1.3p ensures that the limit 
sysyem as £ ^ is not of an Euler type. Also notice that the parameter e in the front of 
the magnetic field H in (|1.2p can be understood as the reciprocal of Alfven number ( [2 7j). 
Under those scalings. system yields 

Pt + div(pu) = 0, 

(pu)t + div(pu (g) u) - ^^Au - AjVdivu + ■^Vp'' = (V x H) x H, (1.5) 
Ht - V X (u X H) = -V X (jy^V X H), divH 0. 

The existence of global weak solutions to (|1.5p has been investigated in Hu-Wang [TB] (and 
in Hu-Wang [T7] for the non-isentropic case). From the mathematical point of view, it is 
reasonable to expect that, as p ^ 1, the first equation in (|1.5p yields the limit: divu = 0, 
which is the incompressible condition of a fluid, and the first two terms in the second 
equation of (jl.Sp become 

Ut + div(u (g) u) = Ut + (u • V)u. 
On the other hand, the incompressible MHD equations read 

Ut + (u • V)u - pAu + Vp = (V X H) X H, 

Ht - V X (u X H) = -V X (i/V X H), (1.6) 
divu = 0, divH = 0. 

Thus, roughly speaking, it is also reasonable to expect from the mathematical point of 
view that weak solutions of p.5p converge in certain suitable functional spaces to the weak 
solutions of (|1.6p as p goes to a constant such as 1 and £ goes to 0, and the hydrostatic 
pressure p in (|1.6p is the "hmit" of {p"^ — l)/e^ in (|1.5p . This paper is devoted to the 
rigorous justification of the convergence of that incompressible limit (i.e., the low Mach 
number limit) for global weak solutions of the compressible isentropic MHD equations. 

In this paper, we shall establish the incompressible limit of p.5p in three types of 
spatial domains: the torus T (in this case, all the functions arc defined on and assumed 
to be periodic with period 27t for all directions, that is, T = [0, 27r]^), the whole space 
M^, and a sufficiently smooth bounded domain G R^, N = 2,3. The study in the 
bounded smooth domain with no-slip boundary condition on the velocity is much harder 
than that in other two cases, because in bounded domains, there are extra difficulties 
arising from the appearance of the boundary layers, and the subtle interactions between 
dissipative effects and wave propagation near the boundary, and hence requires a different 
approach. We remark that the incompressible limits for compressible isentropic Navier- 
Stokes equations have been investigated in [25] for the whole space and the periodic 
domain using the group method, and in [9] for a bounded domain. These results have 
been extended by others, such as [H O [HI [Ml [M] ■ We also notice that in (TS], convergence 
results were proved for well-prepared data as long as the solution of incompressible limit is 
suitably smooth. For the case of non-isentropic flows, see [12j [T3| for some recent studies. 
For other related studies on the incompressible limits of viscous and inviscid flows, see 
[Il[7l[IIl[l6l[l9l[20l[23l[28l[29l[3ll[32] and the references in [12]. Comparing with those works 
on the compressible Navier-Stokes equations, we will encounter extra difficulties in studying 
the compressible MHD equations. More precisely, besides the possible oscillation of the 
density, the appearance of the boundary layer and the interactions between dissipative 
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effects and wave propagation, the appearance of the magnetic field and the coupling effect 
between the hydrodynamic motion and the magnetic field should also been taken into 
considerations with new estimates. We will overcome all these difficulties by using the 
group method, Strichartz's estimate, and the weak convergence method to establish the 
convergence of weak solutions of the compressible isentropic MHD equations (|1.5p to weak 
solutions of the incompressible MHD equations (|1.6p as the density goes to a constant and 
e goes to in the periodic case and the whole space case. More precisely, we will show that, 
for any fixed T > 0, in the periodic case, the incompressible part of the velocity strongly 
converges to a divergence-free vector field in L^([0, T], L^(T)) while the gradient part of 
the velocity converges weakly to zero; and in the whole space case, due to Strichartz's 
estimate, the gradient part of the velocity converges strongly to in L^([0, T], L^q^(M^)), 
while the strong convergence of the incompressible part of the velocity only holds in the 
local sense. However, this method does not apply to the case of bounded domains because 
of subtle interactions between dissipative effects and wave propagation near the boundary. 
Instead, we will use the spectral analysis of the semigroup generated by the dissipative wave 
operator, together with Duhamel's principle. Finally, we remark that the incompressible 
flow also can be derived from the vanishing Debye length type limit of a compressible flows 
with a Poisson damping. We refer the interested readers to |31[5]- 

We organize the rest of the paper as follows. In Section 2, we will give the setting of 
our problem and state our main results. In Section 3, we discuss the convergence of the 
incompressible limit in the periodic case. In Section 4, we will investigate the convergence 
of the incompressible limit in the whole space R^. Finally, in Section 5, wc will study the 
convergence of the incompressible limit in the bounded domain. 



2. Main Results 

In this section, we describe the setting of our problem and state our main results. First, 
we denote by P the orthogonal projection onto incompressible vector fields, i.e. 

V = Pv + Qv, with div(Pw) = 0, cm-l{Qv) = 0, 

for all V € L^. Indeed, in view of results in we know that the operators P and Q 
are linear bounded operators in W^'^ for all s > and 1 < p < oo in the whole space or 
bounded domains with smooth boundaries. Second, let us explain the notation of weak 
solutions to the incompressible MHD equations as follows: Given the initial conditions 
Uq G L^, Ho € such that divuo = and divHo = 0, (u, H) is a weak solution of (|1.6p 
satisfying 

u|t=o = Uo, H|t=o = Ho, (2.1) 

where 

u e C([0, T];Ll,,,) n L2([0, T];H'{n)), H e C([0, T]; L^,,) n L'{[0, T];H\n)), 
if for aU T < oo, V e C^{^) with divt/j = 0, and ip G C^{[0, T)), we have 



^(0) 



/ UQipdx + / "(p'it) / u ■ ipdxdt + / 'tp{t) / {uiditpjUj — fi\/u : Vip) dxdt 
Jn Jo Jn Jo Jn 

[ [ V^(V X H) X H • ipdxdt, 
Jo Jn 



lo Jn 

and 

ft r rt 



ipiO) [ Uoifdx + [ ?/''(*) / H • ipdxdt + [ i'{t) [ (u X H) • (V X ip)dxdt 
Jn Jo Jn Jo Jn 

= v I / (V X H) • (V X ip)dxdt. 
Jo Jn 
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For more details as to the existence and regularity of weak solutions to the incompressible 
MHD equations, we refer the readers to [inilSS]- Now, we can state our main results case 
by case. 

2.1. The periodic case. Let us begin with the periodic case. We consider a sequence 
of global weak solutions {p^,u^,ili.) of the compressible MHD equations (|1.5|) in T and 
assume that 

Pe e L°°([0,T];L^(T)), u, e L^i[0,T]; H\T)), 
Peluep e L^{[0,T];L\T)), p,u, e C (^[0, T]; L^J^, 

H, eL2([0,r];i7i(T))nC([0,r];iLafe(T)), 
for all T S (0,oo), where C{[0,T]; L^^^jS) denotes the functions which are continuous with 
respect to t S [0, T] with values in LP endowed with the weak topology. We require (|1.5p 
to hold in the sense of distributions. Finally, we prescribe initial conditions 

Pe\t=o^p"e, p.u,|t=o=™?-P°uO, H,|t^o = HO, (2.2) 
where p° > 0, p° e i'^(T), e i27/(7+i)(T)^ mO = on {p° = 0}, p°|u^P £ L^T), and 
€ L'^(T). Furthermore, we assume that y^u^ and H° converge weakly in to Uq 
and Hq respectively, and that we have 

' ^ (p°iuOp + |H°p) dx + -^-^ / {{p^r - jpUpir-' + (7 - i)(prr) < 



2 7t'"^' e2(7-l) 
= (27r)"^ / p°d.T ^1, as £ ^ 0, 



(2.3) 



where and hereafter C denotes a generic positive constant independent of e. Notice that 
(|2.3p implies that, roughly speaking, p° is of order ]o7 + 0{e). We assume finally that the 
total energy is conserved in the sense: 



where 



E^{t) + [D^{s)ds < E°, a.e t € [0,T], 
Jo 



(2.4) 



dx. 



and 



1 

where $7 is equal to T in the periodic case, and later is the whole space or a bounded 
domain. 

We now recall the results in [18] which yield the existence of such a solution with the 
above properties precisely as 7 > ^, for = 2, 3. We state the following theorem: 

Theorem 2.1 (The periodic case). Assume that {(p^, u^, HE)}e>o is a sequence of 
weak solutions to the compressible MHD equations (|1.5p in the periodic domain T with 
initial data {(p^, u^, Hg)}£>o, satisfying the conditions (j2.2[) - (|2.4[) and •y > , N ~ 2,3. 
Also assume that (u, H) e [L'^{[0,T]; H^{T))r\ L°° {[Q,T]; L'^{T))]'^ is a weak solution to the 
incompressible MHD equations (|1.6[) with initial data u\t=o = Puq andU\t=o ~ Uq. Then, 
for any finite number T, up to a subsequence, the global weak solutions {(pe, u^, He)}e>o 
converge to (u,H). More precisely, as e —^ 0, 

Pe converges to 1 in C{[0,T]; L~' [Vl)); 
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Pue converges strongly to u in L {[0,T]; L^{T)), for all 1 < p < 



7V-2' 

Qug converges weakly to in L'^{[0,T]; H^{T)); 
converges to H strongly m L'^{[Q,T]; L'^{T)) and weakly in L'^{[0,T]; H^{T)), 
where, for convenience we will denote oo by if N = 2 in this paper. 

2.2. The whole space case. Next, we turn to the whole space case. For the convenience 
of presentation, we only discuss the case when a = 1. In order to define weak solutions in 
the whole space, the following special type of Orlicz spaces i^(il) are needed (see Appendix 
A in [24]): 

LP{n) = {/ e Ll,{n) : fxUfKn} ^ L'^i^), fX{\f\>,} e L^in), for some 77 > 0} , 

where x denotes the characteristic function of a set. We consider a sequence of weak solu- 
tions {{pe, Ue, He)}e>o in the whole space with initial data u°, H°)}e>o, satisfying 
the same conditions (|2.2p and (|2.4p as in the periodic case. In addition, the weak solutions 
{(pe, Ug, He)}e>o satisfy the following conditions at infinity: 

-^1, Ue 0, He 0, as |x| — * 00, 

and 

U (p^|uOp + |H^p)d. + -j^^ / ((p^r-7P° + (7-l))rf-<C. (2.5) 

As pointed out in [TH], one can show that for any fixed e > 0, there exists a global weak 
solution (/9e,Ue,He) to the compressible MHD equations (|1.5p defined by 

Pe-l€L-([0,r];i^(R^)), 

Vp7ue e L°°{[0,T];L\R^)), 

Vue e L\[0,T];L\m.^)), 

He e L2([o,r];ij-i(K^'))nL°°([o,r];i2(R^)), 

satisfying, in addition, 

PeUeeC([0,T];L^;/(^+i)(R~)), 

PeGC([0,T];LUR^)), 
if 1 < p < 7 for all finite number T. 

Now we are ready to state our result in the whole space as follows. 

Theorem 2.2 (The whole space case). Assume that {(pe, Ue, He)}e>o is a sequence 
of weak solutions to the compressible MHD equations p.Sp in the whole space with 
the initial data {(p^, u^, H^)}e>o> satisfying the conditions (|2.2p . (12. 4p . (|2.5p and j > ^! 
TV = 2,3. Also assume that (u,H) G [L^{[0,T]; H^R^)) n L°°{[0,T]- L'^iR^))]^ is a weak 
solution to the incompressible MHD equations p.6p with initial data u\t=o = Pu.a and 
H|(=o = Hq. Then, for any finite number T, up to a subsequence, the global weak solutions 
{(pe, Ue, He)}e>o Converge to (u, H). More precisely, as e 0, 

Pe converges to 1 in C{^,T]; L~' {^l)); 

2na ^1. TP mNv ^ - 2Ar 



PUe converges strongly to u in L ([0, T]; L^^^(R )), for all 1 < p < — -, 

27V 

Que converges strongly to in i^([0, T]; L''(R^)), for all 2 < q < 



TV - 2' 

He converges to H strongly in L^{[Q,T]; L'^{R^)) and weakly in L'^{[0,T]; H^{R^)). 
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2.3. The bounded domain case. The third case we will address in this paper is the 
incompressible limit in a bounded domain f2. For the convenience of presentation, we also 
only discuss the situation when a = 1. In order to state precisely our main theorem, 
we first introduce a geometrical condition on H (cf. [9]). Let us consider the following 
over-determined problem 

— Ai/j = Xijj in Vl, — — = on 3f2, and -0 is constant on (2.6) 

an 

A solution to (|2.6p is said to be trivial if A = and V' is a constant. We say that Q. satisfies 
the assumption (A) if all the solutions to (|2.6p are trivial. In the two dimensional space, 
it is proved that every bounded, simply connected open set Q, with Lipschitz boundary 
satisfies (A). 

We consider a sequence of weak solutions {{pe, Ug, H£)}e>o in a bounded domain 17 with 
initial data {(pg, u°, H°)}e>o and boundary condition 

uelan = 0, Helan = 0, (2.7) 
satisfying the same conditions (|2.2p and (|2.4p as in the periodic case. And the initial data 
of the weak solutions {(pe, u^, He)}j>o satisfy 

\ {pMf + dx + ^^^^ J^^ {{p'^r - 7P° + (7 - 1)) dx < a (2.8) 

As shown in [T^, for any fixed e > 0, there exists a global weak solution (p^jU^jH^) to 
the compressible MHD equations (jl.Sp defined by 

PeeL-([o,r];i^(f7)), 

V^u, eL°°([0,r];L2(f7)), 
Vu,€L\[0,T];H\n)), 

satisfying, in addition, 

PeUeGC([0,T];L27/(7+l)(f^))^ 

PeeC([o,r];iL(f^)), 

if 1 < p < 7 for all finite number T. 

Our main result in bounded domains reads as follows. 

Theorem 2.3 (The bounded domain case). Assume that {{p^,u^,Hg)}^yQ is a se- 
quence of weak solutions to the compressible MHD equations (|1.5p in a bounded domain 
with initial data {(p^, u°, H^)}£>o and boundary condition (j2.7p . satisfying the conditions 
^M, (El), (EH) and > ^, N = 2,3. Also assume that (u,H) e [L^{[0,T]; H^{n)) n 
L°°([0, r]; i^(r2))]^ is a weak solution to the incompressible MHD equations (jl.6p with ini- 
tial data u|t=o P^o o,nd H|t=o = Hq and boundary conditions \i\dn — and H|gsi = 0. 
Then for any finite number T, as e goes to 0, the global weak solutions {(pe, u^, He)}e>o 
converges to (u, H). More precisely, as £ — s- 0, 

Pe converges to 1 in C([0, T]; L'''(il)); 
Uj converges to u weakly in L^{Q x (0,T)) and strongly if Q satisfies (A); 
He converges to H strongly in L'^{[0,T]; L'^{fl)) and weakly in L'^{[0,T]; H^{il)). 

Remark 2.1. In fact, we will split the eigenvectors {^'fe,o}fceN of the Laplace equation with 
Neumann boundary condition into two classes: those which are not constant on dfl will 
generate boundary layer and will be quickly damped, thus converge strongly to 0; those 
which are constant on for which no boundary layer forms, will remain oscillating 
forever, and lead to only weak convergence. Hence, if (A) is not satisfied, u^ will in 
general only converge weakly and not strongly to u. In particular, in the bounded, simply 



LOW MACH NUMBER LIMIT 



7 



connected open set J7 C with Lipschitz boundary, the boundary layer will always be 
generated, and hence will strongly converge to zero. 



3. The Periodic Case 

In this section, wc will prove Theorem 12. II 

3.1. A priori bounds and consequences. We first deduce from (|2.4p and from the 
conservation of mass that we have for almost alH > 0, 



< 



+ [ [ (Mel^'Uep + Ae(divue)^ + J^£|V X Hep) dxds 
Jq Jt 

I (^p°|uOp + IHOp + -^^-^ ((pO)- - jpUpiy-' + (7 - dx < C. 



(3.1) 

From this inequality we see that pejugp, IHgp and i^pj — 'yp^{'p^)'^^^ + (7 — l){'p^)'^) are 
bounded in L^{[Q,T]; L^{T)) and that Du. and V x are bounded in L'^{[0,T]; L'^{T)). 
In particular, we see that pe is bounded in L°°{[0,T]; L^{T)) for all T S (0,00) due to the 
fact that for e small enough, p7 S (l: §) ^^'^ thus for all 6 > 0, there exists some 77 > 
such that 

x^ + i^-l){p;y --,x{p;y-^ >r]\x~p:\'' if \x--p;\>5, x>0. (3.2) 

As in [25] u^ is bounded in L'^{[0,T]; H\T)) for all T G {Q,oo). In fact, we deduce from 
Holder and Poincarc's inequalities that we have for all T e (0, 00) 



JT 



{2n) 



-N 



u^dx 



dxdt < C||pj||ioc([o,T]:L-')II^U£|lL2([0,T]:L2) < C, 



hence, in view of the above bound on pejugp, we get 
rT ^ / !• \ 2 



C> 



Pe 



Wpdx dxdt 



p'idx 



Urdx dt 



Since ()2.3p implies that converges to 1 in measure, and hence, up to a subsequence, in 
L^(T), thus, wc can deduce a bound on u^ in L^([0,r];L^) by using Poincare inequality 
again. Indeed, we have 



\\iA^dxdt < 2 



-N 



u^dx 



dxdt + 2{2tt)' 



dxdt 



< C 1 



JT 



WuJ^dxdt < C. 



From now on, we assume that, up to a subsequence, Ug converges weakly to some u in 
L^{[0,T];H\T)) for all T > 0. On the other hand, the bound on in L°°{[0,T]; L^{T)) 
and the bound on VH^ in L^([0, T]; L^(T)), combining the following Gagliardo-Nirenberg 
inequality 

i 1 

"""l§([0,T];L4(T)) - II"IIl-([0,T];L2(T))II^"IIl2([0,T];L2(T))' 

imply that is bounded in L^/^{[0, T]; i'^(T)), and also we can assume that converges 
weakly to some H in i^([0, T]; (T)) with divH = 0. Finally, from the induction equation 
in (|1.5p . we sec that 

dtHe = V X (Ue X He) - V X (l/^V X H,) 
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is bounded in L^/'^([0, T]; iJ-i(T)), because the fact that is bounded in L'^{[0,T]; L'^{T)) 
implies that x and i^^V x are bounded in L*/'^([0, T]; ^^(T)). Then the Aubin- 
Lions compactness Lemma (see |24| ) implies that 

He^H, strongly in L^^'^ {[0,T]; L^{T)). 

Moreover, this, combing with the uniform bound on Hg in L°°([0, T]; L^(T)), implies that 
He converges strongly to H in L^([0, T]; L^(T)). Therefore, by a standard argument, we 
deduce that the limits u and H satisfy the induction equation in the sense of distributions, 
and also the nonlinear term (V x H^) x in the second equation of (jl.Sp converges to 
(V X H) X H in the sense of distributions. 

Next, we claim that converges to 1 in C{[0,T]; L^{T)). Indeed, in view of (j3.1|l and 
(13.21). we have 



sup / \pe - l^dx < S'^{2tt)'^ + Csup ( / X{\p,-l\>5}\Pe - Pel'^dx 

t>o Jt t>o \Jt 



■C\p-e 



(-12 

< (27r)^J^ + — +qp7-ir, 
V 

and we conclude the claim upon letting first e go to and then J go to 0. 

Now, we show from the previous bounds that divu,. converges weakly to in L'^{[0, T];L'^{T)) 
and that Pu^ converges to u = Pu strongly in i^([0, T]; L^(T)), and thus by Sobolev 
imbedding in L'^{[Q,T]; L"^) for all 2 < q < These facts imply that converges 

weakly to in L^([0, T]; iJ^(T)). Indeed, since converges to 1 in C{{0,oo);L'^{T)) and 
7 > we deduce from pTS]) that divu. converges weakly to in L^([0, T]; L^(T)). The 
second part is proven by observing first that we project (jl.Sp onto divergence-free vector- 
fields: 

dtP{peU,) + P[div{p,u, (E) u,)] - p.APu, = P((V X H,) x H,). (3.3) 

Noticing the fact that the operator P is bounded in all Sobolev space W'^'P for all s G 
[0, oo) and 1 < p < oo and the preceding bounds, p.3p yields a bound on dtP{p^\ie) in 
L\%T]-H-\T)) + L^{[0,T]-L\T)) + L\[Q,nH-HT)), hence, in L\[0,T]- H-\T)). 
In addition, P{peVL^) is bounded in L°= ([0, T]; (T)) n L2([0, T]; L''(T)) with 

1 1 N -2 



r 7 2N 

Next, we will need the following compactness Lemma (cf. Lemma 5.1 in [H]): 

Lemma 3.1. Let gn, hn converge weakly to g, h respectively in L^^ (0, T\ L^^), L'^ (0, T; L''^) 
where 1 < pi,P2 < oo, 

1 11 1 _ _^ 

Pi qi P2 92 

Assume in addition that 

is bounded in (0, T; for some m>0 independent of n, 

dt 

and 

\\hn (• +f'^)llL5i(o,T;L92) ^ as 1^1 ^ 0, uniformly in n. 
Then converges to gh in the sense of distributions in x (0, T). 

Applying this lemma with the previous bounds, we deduce that P(peUg) ■ Pu^ con- 
verges in the sense of distributions to |up. We then conclude easily that Pu^ converges 
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in L2([0, T]; L2(t)) to u upon using the weak convergence of Pu^ to u in L'^{[0, T];L'^{T)) 



and remarking that we have 



f f (|Pu,|2 - P{p,u,) ■ Pu,) dxdt 
Jo Jt 



< C\\pe 



|C([0,T];L^) 



IL2([0,T];L = )' 



with s = < since 7 > ^• 



We conclude this first step by showing the foUowing bounds vahd for aU i? e (l,oo) 

ll<y5e||L°°([0,T];L2(T)) < C* if 7 > 2, 



IIV'eX{pe<«}IU==([0.T];L2(T)) < C if 7 < 2, 



(3.4) 



if 7<2, 



\WeX{p^>R]\\L^([a,T]-L'i(T)) < Cei 

where we denote the density fluctuation by (/j^ = ^ (p^ — p7) • These bounds are deduced 
immediately from the following straightforward inequalities: for some v > Q and for all 
X > 0, 

a;'' - 1 - 7(x- 1) > - 1|2 if 7 > 2, 

a;'' - 1 - 7(x- 1) > - 1|2 if 7<2 and x < R, (3.5) 
x'^ - I - -l{x - I) > v\x - if 7<2 and x>R. 

3.2. The weak convergence of Qn. The proof of the weak convergence of Qu is similar 
to that in Lions-Masmoudi [25] , thus we only briefly describe the main idea from ^5] . We 
provide here first a formal proof of the passage to the limit, next the main difficulty, and 
finally the strategy of proof used in order to circumvent that difficulty. 

We thus begin by an informal proof. It is not difficult to check that the main difhculty 
with the passage to the limit lies with the term div(/OeUg ® u^) and more precisely with 
the term div(p£(5(ue) ® Qu^) since the strong convergence of Pu^. Formally, this term 
should not create an obstruction since in view of the continuity equation in we can 

rewrite the term 9t(peUe) + div(p£Ue ® Ug) as PedtUg + /Oe(ue • V)u£:, which corresponds 
to the term dtu + (u • V)u in incompressible MHD equations (|1.6p . Next, the dangerous 
term [(Qu^) • VjQu^ can be incorporated in the pressure p at the limit since Qu^ = V-0e 
for some ip^, and then 



[(Que) • V]gu, = V 



Next, we need to write down rigorously the proof of the convergence. First, we introduce 
the following group {£.{t),t G R} defined by e*^ where L is the operator defined on Vq x 
{VT, where P^, = {0 € P', / = 0}, by: 



/ divu\ 



for 6 > 0. 



(3.6) 



We remark that e*^ is an isometry on each H'^ x (iJ*)^ for all s e i? and for all t, endowed 
with the norm ||((/),u)|j = (||0|||j^s + ^or details, we refer the reader to pS] . 

For convenience, in the sequel, we will denote by Ci (£2) the first (the second) component 
of the operator C, respectively. 

We next claim that 'C(— |)(g^J°^ ^) is relatively compact in i^([0, T]; i/^") for some 

n e (0, 1). To this end, we need to prove first that (q(^\, )) is bounded in L^{[0,T]; H~") 
for some s £ (0, 1) and that dt{C (-|) (qJ/u,))} is bounded in L^{[0,T]; H"'') for some 
r > large enough. Our claim then follows from Aubin-Lions compactness lemma by 
choosing n in (s, 1). 

Since we know that cp^ is bounded in L°°{[0,T]; LP) where p = min(2,7), by Sobolev 
imbedding theorems, we know that ip^ is bounded in L^([0, T]; i/^**) for some s € (0,1]. 
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And, we also deduce from the previous subsection that psU^ and thus Q{peUe) is bounded 
in L2([0,r];L«) with 

1 _ 1 N -2 
g ^ 7 ^ 2N ' 

Therefore, C (-|) (^(J/u^)) is bounded in L^{[0,T]; H^") for some s G (0, 1). 

In order to get the uniform bound on dt{C (—7) (g(J u ))}' project the second equa- 
tion of p.Sp into the space of gradient vector-fields and we find 
dtQip^u^) + Q[div(peU£ » u^)] - + Ae)VdivUe4- 

(pj - 7P.(pr)^-' + (7 - + °^^^i''' v(pe - p7) (3.7) 

= g[(VxH,)xH,]. 
Hence, we can write the first equation (jl.Sp and (|3.7p as 

£^ +divQ(p,u,) = 0, 

where 6 = 07(^7)'''"^, and 

=(yUe + Aj)VdivUe - Q[div{p^u^ (X) u^)] 

- aV (^^(pj - ^PeCiPeW-' + (7 - l)(pr)^)] + Q[(V X H,) X H, 

Since b goes to 07 as e goes to 0. we will ignore the dependence of 6 on e hereafter. Now, 
we set 

then we have 



£ \Q{PeU^)J e \Q{peUs 



d_f \ l/divQ(p,u,; 

9t \m^J ^ \ £ J \ dt \Q{p^VL.) J e \ bVipe 



A/0 



where is bounded in L'^{[Q,T]:H-'^{T)) + L'^{[0,T]-W-^-'''^{T)) for all (5 > 0, and 

4 + 1. Thus 



hence, is bounded in L'^{[0,T]; R-'' {T)) for all r > f + 1. Thus is bounded in 



L2([0,T];iJ-'-(T)). 

We deduce from the compactness of {ip^i'iTT-e) that we may assume without loss of gen- 
erality that (i/jejTOe) converges in L^([0, T]; iJ^") to some {ip,m). Since Pm^ = 0, we also 
have Pm = 0. Similarly, J ip = 0. Hence, we have 

^ - c(-] ('^] +re, re^O in L^{[Q,T]; R-") as e ^ 0. (3.8) 



^<9(PeUe)y \£j \my 

Finally, following the argument of Step 4 in Section 3 in [25], one can show that (p, 
m £ L^([0, T]; iv^(T)) and div(p£Ue (E) u^) — div(ue (E) v^) converges to div(u (E) u) in the 
sense of distributions, where = £2(7)(^). Moreover, following the argument of Step 
5 in Section 3 in |25| . we can show that div(i;£ (E v^) converges to a distribution which is 
a gradient. Note that the magnetic field does not affect the argument of convergence of 
div(peUg E) Ug) — div(z)e E) Ve) and div(t;£ (E v^) because the magnetic field H does not affect 
the integrability of based on our estimates, thus we only state those convergence results 
without proof. We refer the reader to [53] for details. 

This finishes the proof of Theorem 12.11 
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4. The Whole Space Case 

In this section, wc prove Theorem 12.21 The idea is taken from [5]. Before we start, we 
introduce homogeneous Sobolev spaces for 1 < p < oo and s € R defined as usual by 

W'-PiR'^) = {-A)-'/^LP{R^) and H'{M.^) = Ty"-2(M^), 

where A is the Laplace operator. 

Let us denote by C € C^(R^) a smoothing kernel such that C ^ 0; /rn (dx — 1, and 
define Caix) = a~'^({x/a). The following estimate will be useful in this section (cf. [5]): 

\\f-f*U\L.<Ca^-"\\Wf\\L2, for all f e H\ (4.1) 

where 

2, and a = N { 



and for 1 < p2 < pi < oo, s > and a e (0, 1), we have 

||.9*ailL.i(K«) <Ca-^-^(i/^'^-i/P^)||5l!M/-...(K«). (4.2) 

4.1. A priori estimates and consequences. Most of the arguments developed in the 

periodic case can be adapted to the whole space case. First, we obtain bounds on Du^ in 
L2 ( [0, T] ; L2 (rW ) ) ^ on V X He in L2 ( [0, T] ; L2 (M^ ) ) and on I Ue P , and (p7 + _ 1 ) _ 

7Pe) in L°°([0, T]; Li(R^)). The bound on in L^{[0, T];L^{R^)) follows from ^ and 
the following observation: 



J^^ (^IPs - l|'X{|p.-l|<l/2} + -^\Pe - irX{|p.-l|>l/2}) < C, 

and thus, in particular. 



(4.3) 



/ \Uefx{p,<l/2}dx 



where 



<C+( f X{p.<i/2}dx) / |u,|27'dz) 
<c{l + (meas(|p, - 1| > l/2))i/7||u,|||» pu,||'^i^- 
<c{l + s'/^u,\\l^Du,\\ir'^' 

..N-2 1 



2 ^ ' 2N 2y 

We then complete the proof of our claim using the bound on Du^ in L^([0, T]; i^(]R^)) 
and the classical Young's inequality. Moreover, if we define the density fluctuation as 

e 

then, it is bounded uniformly in e in L°°([0, T]; ij) with k — min{2,7}. Furthermore, if 
we write 

where 



= UeX{|pe-i|<i/2} and = UeX{|p.-i|>i/2}, 

then, we have 

sup/ lu^pda; < 2sup / pe\us\'^dx < C, 
t>0 JR« t>o Jk" 
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and for p < k when N ~ 2, p ~ 2k/3 if = 3, 

/ \nl\^dx<c[ |p,-irx{|p.-i|>i/2}|u.pda; 



2 



< C\\{pe - l)X{|p,-l|>l/2}|lLc»([o.T]:L«(R"))ll"ellL2»/(«-p) 



hence, by Young's inequahty, is bounded in i°°([0, T]; L^(M^)) and u^e ^ is bounded 
in L'^{[Q,T]-L'^{m.^)), where /3 € (0, 1) if iV = 2 and /3 = 2/3 if = 3. 
RecaUing that 7 > iV/2, we deduce that is bounded in 

L2([0,r];i^(R^)ni27/(7-i)(]RA^)), 

Hence, we have 

IIVeUe||^2Qo,T];LV3(K«)+L2«/('' + i)(R") < C. 

Therefore, using Sobolev's imbedding, we deduce 

ll</'eUe||L2([o.T]:_ff-i(R")) ^ C. 

Finally, we already know that 93° is bounded in L2(M^), hence in H~^{R.^), since 
7 > iV/2. On the other hand, m° can be rewritten as 



"^e = ^VP?X{|p?-l|<l/2} + ^ , 0_li VlPg - l|X{|pO_i|>i/2}. 

This implies that m° is bounded in L2(M^) + L2«/(«+i)(M"), and hence in i7-i(M^). 
Therefore, (^f,) is bounded in ff~^(M^) uniformly in £. 

4.2. Strong convergence of Qu^ to 0. We now prove that the gradient part of the 
velocity Qu^ converges strongly to 0. More precisely, we claim that Qu^ converges strongly 
to in L2([o^t];LP(R")) for all p G (2, . Indeed, let us first observe that the 

compressible MHD equations can be rewritten in terms of the density fluctuation the 
momentum = PeU^ and 0^ ~ (m^) follows 

where the wave operator L is defined on (I?'(R^))^+^ with values in (I?'(R^))^+^ by 

/ div ra\ ( 



and 



MeAul _^ AeVdivul - div(TOe (g) Ue) - j^V(/97 - 1 - 7(pe - 1)) + (V X H^) X H, 





/ieAu2 + AeVdivu^^ 
Using DuhameVs formula, we deduce that 



,it) = C (^^ QcHi + f^C {^-^^ {QFl{s) + QF!{s))ds. 



Here we used the fact that Q and C commute, since Q and L do. 

At this stage, the following Strichartz's estimates from [8] are useful: 
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Lemma 4.1. For all s > 0, we have 



( - ) QVo 



Li{(0,oo)-W' 



<Cei/'||^o||//-MK~)' 



V e 



Qij}{s)ds 



<C(l+T)£l/«||V'|U,([O.Tl;ff- 



L9([0,T];iy- = -<'.p(K")) 

/or aZ/ (p, q) G (2, oo) x (2, oo) and cr G (0, oo) such that 

2 , , /I 1\ iV+ 1 



(4.4) 
(4.5) 

(4.6) 



Now, we choose p G (2, jf^j )' 9 ^ ^^^"^ ''^ '= °°) given by (|4.6p . One can deduce 

that 



Hence, 



IQUel < IQUe - QUe * Cq I + £|Q(Ue(;5e) * Ca I + IQ™e * Ca|- 
l|QUe|l2,2([o^T]:LP(K«)) < C'"^"^*'^^^""^^''' I! VUe||i2([o,T];L2(KW)) 

+ ea-i-~(V2-i/p)|,^^u^||^,^j^_^j^^_,^^„^^ 

+ IIQnT-c * Ca||L2([0,T];LP(R"))- 

From the estimates in previous subsection, we know that is bounded in L°°{[Q, T]; iJ~^") 
for all So > N/2 + 1. On the other hand, we deduce from the uniform bound on u^e^^ in 
i2([0,T];L2) that e-^^F^ is bounded in ^^([O, T]; iJ-2(M^)). Then, using LcmmaHIIl we 
obtain, for all ?/ > small enough, 

WQm^ * Ca||L2([0,T]:LP(R")) 
* ) i'e 

L<!([0,T];W-i-<'-P(R")) 
T 



L<!([0,T];W-'J-«/2-l-''.P(R«)) 



-2-Af(l/2-l/p) 



V e 



L2([0,T];_ff-2) 

< Cra-'-^e^/'^UlWH-^ + CTa-^/2-^-'^-''e^/''||F,i|L.([o,T];i/-.,-»/-i) 

Next, fixing a > and letting e go to zero, wc obtain 

limsup|iQu,|U2([o,T].LP(R«)) < Ca^-''^^'^-^'P\ 

where C is independent of e and a. Noticing that 1 — N{l/2 — 1/p) > 0, wc finally get, by 
letting a ^ 0, 

limSUp||Qu^||i2([o,T]:LP(R«)) = 0. 

This implies that Qu^ strongly converges to in L'^{[0,T]; LP(R^)) for aU 2 < p < 
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4.3. Strong convergences of Pug and Hg. In the previous section, we proved the strong 
convergence of the gradient part of the velocity to 0. In order to complete the proof of The- 
orem [521 we are left to show the convergence of the incompressible part the velocity, Pu^, 
the convergence of the density, and the convergence of the magnetic field. This can be done 
by using the classical compactness arguments in [23 , or equivalently by looking at the 
time-regularity properties of Pu^, see [5]. Indeed, following the argument in the periodic 
case steps by steps, we obtain the strong convergence of to 1 in C([0, T]; Lj'^^(M^)) and 
the weak convergence of Pu^ to u in L^([0,T]; Moreover, we also can show that 

Pug converges to u in L^([0, P]; L^(Ph)) for all R € (0, oo). Here, we denote by Bn the 
open ball centered at of radius R. 

Finally, similarly as in the periodic case, the bound on in L°°([0, P]; P^(R^)) and 
the bound on VH^ in P^([0, P]; P^(M^)), combining Sobolev's inequality and interpolation 
theorem, we know that H. is bounded in P*/3([0, P]; P^(R^)), and also we can assume that 
converges weakly to some H in p2([o^T];iJi(]R^)) with divH = 0. Finally, from the 
induction equation in (jl.Sp , we deduce that dtH^ is bounded in P*/^([0, P]; iJ-i(M^)), due 
to the fact that is bounded in P^([0, P]; P^(R^)). This property, combining Aubin-Lions 
compactness Lemma, implies that converges strongly to H in L®/^([0, P]; P^q^(R^)). 
Moreover, the uniform bound on in P°°([0, P]; i^(R^)) implies that converges 
strongly to H in L'^{[Q,T]; Lf^^{R'^)). Therefore, by a standard argument, we deduce that 
the limits u and H satisfy the induction equation in (jl.ip in the sense of distributions, 
and also the nonlinear term (V x H^) x in the second equation of (jl.Sp converges to 
(V X H) X H in the sense of distributions. For a detailed statement of the above argument, 
we refer it to the argument surrounding the convergence of the magnetic field in section 
3.1. 

The proof of Theorem 12.21 is complete. 



5. The Bounded Domain Case 



In this section, we will prove Theorem 12.31 by the spectral analysis of the semigroup 
generated by the dissipativc wave operator. Before we start, we introduce the eigenvalues 
{^k olfceN (-^fc,o > 0) and the eigenvectors {^fc,o}fcGN hi L'^{fl) with zero mean value of the 
Laplace operator satisfying homogeneous Neumann boundary conditions: 

-A*fe,o = A^o*fe,o inf^, ^^ = on dQ. 

on 

Notice that, by Gram-Schmidt orthogonalization method, it is possible to assume that 
{^fc,o}feGN is a orthonormal basis of L^(ri) and that up to a slight modification, if Afe,o = Xi.o 
and k I, then 

/ V*fc,o • V*,,od.s = 0. 
Jdn 

Next, we recall from the previous section that we can deduce similarly that 

sup \\pe - l||LT(n) < Ce'^/'' and sup \\p^ - l\\L-{n) < Ce, 

t>0 t>0 

where k ~ min{2,7}. And similarly to the whole space case, we will split 

Ue = + U^, with ul = U,X\p,-l\<l/2, = UeX|p,-l|>l/2, 

which satisfy 

sup / |u^pd2;<2sup / pe\ue\'^dx < C; 
t>o Jn t>o JQ, 



and 



lUelli^(n) < 2 / \Pe~ l\\ne\^dx < Ce|j U, || (j,, < Ce|VUe|i.(o). 
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Therefore, is bounded in L°°{[0, T];L^{n)) whereas u^e-^/z ig bounded in L'^{nx{0, T)), 
and hence, is bounded in L^(r2 x (0,r)). Also, in this section, we denote the density 
fluctuation by 

and the momentum by ~ PeUg. 

5.1. Strong convergence of Pu^ and H^. Following the argument in the periodic case 
step by step, up to the extraction of a subsequence, we then obtain the strong convergence 
of pe to 1 in C([0, T]; LT(17)), the strong convergence of Pu^ to u = Pu in L'^{n x (0, T)), 
and the weak convergence of Qu^ to in L'^{[Q, T]; H^{^l)). Thus, the continuity equation 
in (jl.Sp holds in the sense of distributions. 

Similarly to the periodic case, the bound on in L°° {[0,T]; L'^{n)) and the bound on 
VHg in L^([0, r]; i^(r2)), combining Sobolev's inequality and interpolation theorem, we 
know that is bounded in L^/^{[0, T];L^{n)), and also we can assume that converges 
weakly to some H in L'^{[Q,T]; H^{n)) with divH = 0. Also, from the induction equation 
in p3)) . we deduce that dtU^ is bounded in L**/'^([0, T]; iJ-^(O)), due to the fact that 
is bounded in L^([0, T]; This property, combining Aubin-Lions compactness 

Lemma, implies that converges strongly to H in i*/^([0, T]; L^(r2)). Moreover, the 
uniform bound on in L°°([0, T]; L^(Jl)) implies that converges strongly to H in 
i^([0, T]; L^(il)). Therefore, by a standard argument, we deduce that the limits u and H 
satisfy the induction equation in (jl.ip in the sense of distributions, and also the nonlinear 
term (V x H^) x in the second equation of (|1.5p converges to (V x H) x H in the sense 
of distributions. Therefore, in order to prove Theorem 12.31 it only remains to study the 
convergence of the gradient part of the velocity Qu^. 

5.2. The convergence of Qu^. The argument for the convergence of Qu^ in this subsec- 
tion follows the lines in [S], except the argument for the magnetic field. For the reader's 
convenience and the completeness of the argument, we provide the details here. For this 
purpose, first, we discuss the spectral problem associated with the viscous wave operator 
!/£ in terms of eigenvalues and eigenvectors of the inviscid wave operator L, where the wave 
operator L and are defined on 'D'{n) x I?'(ri)^ by 

, ■0 \ f divm^ 

and 



m/ V V* 



my Vm/ \/ieAm + AjVdivm 
The eigenvalues and eigenvectors of L read as follows 



fe,0 



fe,0 - I _ _i_ V^k o 

= ±iAfc,o0^_o ^' "^fc,o • n = on dQ. 
In the following steps, the following information on the approximating eigenvalues and 
eigenvectors for the operator is crucial: 

Lemma 5.1. Let O be a bounded domain in and let k > 1, M > 0. Then, there 
exists approximate eigenvalues iA^^j^j and eigenvectors (t^keM ~ \ \^'") oj such that 

with 

i\k.e,M = ±i^fc,o + iA^iV^+ 0(e), where 7^e(iAfe J < 0, 
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and for all 1 < p < oo, we have 

Proof. For the construction in detail, wc refer the readers to [5]. □ 
Remark 5.1. Due to the construction in [9], indeed, we have 

l±i 



' Me 
2^^1,0 -Ida 



fe,0| 



■ ds. 



Remark 5.2. We notice that the first order term iA^j^ clearly yields an instantaneous 
damping of the acoustic waves, as soon as 7^e(iAJl) < 0. For this reason, we define / C N 
to be the set of eigenvectors ^k,o of the Laplace operator such that TZe{iX^^) < and 
J = N — /. Observe that when S J, we have ^ =0. For those indices, m^g identically 
vanishes on dQ and therefore satisfies not only • n = but also m^Q = on dQ, hence 
no significant boundary layer is created, and there is no enhanced dissipation of energy in 
these layers. 

Now, we can express Qu^ in the terms of the orthonormal basis | ^Xk'^o" } '^^ L^(ri) 



as 



feeN 



fc,0 



Afc.o 



where the notation (•, •) stands for 



if{x),g{x)) = / f{x)g{x)dx 



We can split Qu^ into two parts Qiu^ and Q2Ue, defined by 



fc,0 



A 



fc,0 



fc,0 



Afc ,0 



and Q2Ue = ^ 



A 



k.O 



'k,0 



Afc,o 



kei ^ ' '"'^ fee,/ 

which respectively correspond to damped terms and nondamped terms. We will prove on 
one hand that Qiu^ converges strongly to in x (0, T)), and on the other hand that 

curldiv((52TOe ® Q^^e) converges to in the sense of distributions, if J 7^ 0, which is equiv- 
alent to say that div(Q2TO£ ® Q-iy^e) converges to a gradient in the sense of distributions. 

Let us observe that in view of the bound on in i^([0, T]; i^Q (fi)), the problem reduces 
to a finite number of terms. Indeed, we have 



fe>M 



fe,0 



dt < 



C 



X2 



I^Ueli2(nx(0,T))J 



i = 1 or 2. 



Hence, recalling that A 



M 



00 as M 



00, we only have to prove that (QiU^, mj^) 



converges strongly to in L^{0,T) for any fixed k, and study the interaction of a finite 
number of terms in div((52Ue ® Q-zU^). On the other hand, we notice that 



and 



e|(Q(^eUe),V*fc,o)| 



< e\\(pe\\L-'in)\\Ue\\^-^ 



fe,0||L°°(O)i 



which goes to in L^(0, T) since 7 > N/2. Hence, we are led to study (Qnie, hi^q). 
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Denote 



we observe that: 



2(Qm„m±o) = /3±^-/3,T, 



so that it suffices to consider the convergence properties of in L^(0,T). Also, we see 
from Lemma [5JJ witli M = 2 that: 



0^0 - -^^,2)1 < Ce"/' (ll^e||L~([0,T];Lnn)) + l|me|| 



L°°([0,T];LWT(q)) 



where 



a = min^l-i 



hence we only have to prove that b"^ ^{t) = {4>e{t)^ 0fc e 2) converges strongly to in L^(0, T) 
when k G I, and study its oscillations when k £ J. 
Notice that 4>£{t) = ((/3£,me) solves 



(5.1) 



where L* denotes the adjoint of with respect to (•, •), and 

"(p,r-7Pe + (7-l) 



5e = — div(me (8) u^) — V 
Taking the scalar product of (|5.ip with 2' obtain 

'^here c^ J<) = (g., m^^^ 3) + '^"H'?^^' -^^,2)- 



(V X H,) X H, 



(5.2) 



5.2.1. The case k e /. From (j5.2p . by Duhamel's principle, we deduce that 

6± (i) = 6± (0) CXp'^t,2*/e ^ ^± CXp'^t,2(*--)A 



(5.3) 







The first term in (|5.3p is estimated as follows 

< C 6^,(0) exp^^('^^i)*/^ 



6± (O)exp 



i^?,.,2*/T 



L2(0,T) 



L2(o,T) 



In order to estimate the remaining term in (|5.3p . we will use the following estimate: for 

< 

^E,2(*-«)A, 



any 1 < p, q < 00 with ^ + i = 1 , we have 



exp 



a(s)ds 



< / exp^"('^'=,i)(*-^)/^|a(s)|ds < C||a|U,(o,T)e*- (5.4) 



We now write |c, | < ci + C2 + C3 + C4, where 



ci(i) = 



C2(i) = 



/ (me (X) Uj)(<) • Vm^^ ^dx 
Jn ' ' 

{p,r - + (7 - 1) 



(t)divm^ 2dx 



C3(t)=£-l|(0e,i?± )l, 



C4(t) 



/ (V X He) X He • Vm^^ 2'^2; 
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Observing that rtie = SipeU^ + u^, we have: 

ci(t) < ||m±^ 2lU~(n)l|u^ + u^llL2(j2)|lVue|lL2(n) 

+ £||</'e||L<»([0,T];L''(f2))||(Mc)^!lL'=/('<-i)(n)l|VmJ^^2||ioo(o) 
< C||ul|U,^([o,T];L2(o))|lVu,|U2(o) +C£l/2||Vu,||2,(j,) +C£l/2||Vu,jU2(o). 

The second term C2 is estimated as 

C2(i) < C'||pe||Loo([o,T];Li(n))(l!*L,2lU°°(iJ) + II ^fc,e,2 II (f^) ) ^ ^■ 

Also, we estimate C3 by 

C3(i) < Jl|i?L,2llL''/(«-)(O)ll</'e||L-([0,T];L~(n)) < Ce^/^-l/^^. 

Finally, we can estimate C4 by 

C4,{t) < ||Vm^^^^2lU~(O)l|He||L-([0,T];L2(o))||VHe|lL2(o) < C|| VH^ || L2(n) . 

Therefore, using the estimate (|5.4p repeatedly, we can conclude that ^ converges strongly 
to in L2(0,T). 

5.2.2. The case k G J. From ((O)) and the fact that AJ^^ = 0, we see that exp±'^'=.o*/^ 
is bounded in L^{0,T) and that its time derivative is bounded in y/eL^{Q,T) + LP{Q,T) 
for some p > 1. It follows that up to a subsequence, it converges strongly in L^{0,T) to 
some element 

Next, since Peix, t) converges to 1 in C([0, T]; L'^{Vl)) and b^ ^{t) are uniformly bounded 
in L^{[0, T]), we deduce that 

in the sense of distributions. Hence, we only need to consider the terms 

biAe^^^^ foraUfc,?e J. 

On the other hand, due to the strong convergence of exp*''^'= °*/'^ in i^([0,T]) when 
fc S J, we can deduce that 

6± 6± - 



= exp'(^^.«-^'.o)*/- exp-'^-«*/^ bk.o exp'^-*/^ 6,,o^^ ® ^ 
- expH^^.o-.o)*/e ,±_(t)6^_(0^ 

at least in the sense of distributions. Thus, we are only left to study the interaction of 
terms 

expi(>.o-A.o)*/. btJt)btJt)^ « (5.5) 

We will finish the analysis of the interaction by two cases. The first case is A^.o = Xi^. 
In this case, the term (|5.5p is reduced to 

whose divergence is clearly a gradient in the sense of distributions, due to the fact that as 
long as Afc.o = A/.o, we have 

div(V4'fc,o ® V^-z^o + V^-z^o ® V^-fc.o) = -A^.oV(1'fe,o*;,o) + V(V*fc,o • V^-z^o)- 
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For the second case, we have Afc,o 7^ A/.o- Under this situation, due to the fact that 
btscit) G L\[0,T]) askeJ, we know that fcL.cW&foscW S LH[Q,T]) for aU kj e J. 
Then, by Riemann-Lebesgue Lemma, we conclude that 

exp'(^'^-^-^'.»)*/e 6t.e(0^Lc(0di - 0, as e ^ 0, 

which imphes that (|5.5p converges to in the sense of distributions. Hence, the finite sum, 
as fc, / < M, 



div PebkA.e^ -\ 



converges to a gradient in the sense of distributions. And hence, 

div(peQ2Ue (g) Q2Ue) 

converges to a gradient in the sense of distributions. 
This completes our proof of Theorem [ 
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